Longitudinal interlayer magnetoresistance in quasi-2D metals 
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The longitudinal interlayer magnetoresistance R zz (B z ) is calculated in strongly anisotropic lay- 
ered metals, when the interlayer band width 4t z is less than the Landau level separation hw c . The 
impurity scattering has much stronger effect in this regime than in 3D metals and leads to a linear 
longitudinal interlayer magnetoresistance R zz oc B z in the interval fkj c > it z >> yTofcl changing 

1/2 

to a square-root dependence R zz oc B z at higher field or smaller t z . The crossover field allows to 
estimate the interlayer transfer integral as t z ~ y/Tohuj c . Longitudinal interlayer magnetoresistance, 
being robust to the increase of temperature or long-range disorder, is easy for measurements and 
provides a useful tool to investigate the electronic structure of quasi-two-dimensional compounds. 

PACS numbers: 72.15.Gd,73.43.Qt,74.70.Kn,74.72.-h 



The investigation of angular and field dependence of mag- 
netoresistance (MR) provides a powerful tool of study- 
ing the electronic properties of various metals, including 
strongly anisotropic layered compounds, such as organic 
metals (see, e.g., Refs. [l|-[4j for reviews), cuprate and_ iron- 



based high-temperature superconductors, (see, e.g., (H-14|) 
heterostructures [lij etc. In layered quasi-two-dimensional 
(Q2D) metals with at least monoclinic crystal symmetry 
the electron dispersion in the tight-binding approximation 
is given by 



e 3D (k) rj e 2 D ~ 2t z cos(k z d), 



(1) 



where the 2D electron dispersion in magnetic field perpen- 
dicular to conducting layers is quantized in Landau levels: 



(2D = £215 (n) = hio c (n + 7) . 



(2) 



where the Landau level (LL) separation huj c — HeB/m*c, m* 
is the effective electron mass, n is the LL number, 7 k, 1/2, 
k z is out-of-plane electron momentum, and d is the inter- 
layer spacing. If the interlayer transfer integral t z ^> Hui c , 
the standard 3D theory of galvanomagnetic properties [l7r- 
flit can be applied, which is valid in the lowest order of the 
parameter Huj c /t z . This theory predicts some special fea- 
tures of MR in Q2D metals: the angular magnetoresistance 
oscillations (AMRO) [20H22] and the beats of the amplitude 
■of magnetic quantum oscillations (MQO.fljj 

In more anisotropic Q2D metals, when t z > hio c , several 
new features appear, such as slow MR oscillations 23|, \24 \ 
and the phase shift of MQO beats between transport and 
thermodynamic quantities. |24l [25j These two effects are not 
described by the standard 3D theory 17 -lj[ because they 
appear in the higher orders in the parameter fku c /t z . The 
monotonic part of MR also changes when fku c /t z ~ 1. 
According to the standard theory, [17J external magnetic 
field along the electric current leads only to MQO but 
does not influence the monotonic (background) part of this 
current. However, the monotonic growth of interlayer MR 
R zz with the increase of longitudinal magnetic field B z was 
observed in various compounds as a general feature of Q2D 
metals. 23ll27H34j Its theoretical description is the aim of the 
present paper. 

In very anisoropic and dirty compounds with t z < To, huj c , 
where Tq = ft/2ro and tq is electron mean free time in 



the absence of magnetic field, this monotonic growth of 
Rzz (B z ) was attributed to the new "strongly incoherent" 
mechanisms of interlayer electron transport, such as metal- 
insulator transition and variable-range electron hopping with 
exponential temperature and field dependence of -Rz Z ,[35[ or 
the interlayer hopping via local crystal defects with in-series 
metallic intralayer transport. [32] The in-plane component of 
magnetic field in the "strongly incoherent" regime leads to 
weaker MR compared to coherent 3D theory. 32, [3^, 37 1 
However, in the most experiments the longitudinal inter- 
layer magnetoresistance was observed together with the pro- 
nounced AMRO and metallic-type temperature dependence, 
which is inconsistent with " strongly incoherent" mechanisms 
of interlayer electron transport, which do not conserve the 
in-plane electron momentum during interlayer electron hop- 
ping. Recently it was shown. [38l44fj | that at very weak in- 
terlayer coupling, hcu c ^> Tq ^> t z , the longitudinal inter- 
layer magnetoresistance has a square-root monotonic growth 
R zz oc V ' B z even within the coherent-tunnelling model. The 
angular dependence of MR also changes in this limit, [HI 
which contradicts the previous common opinion 41| that in 
the "weakly incoherent" regime, i.e. at Tq > t Z) the inter- 
layer magnetoresistance does not differ from the coherent al- 
most 3D limit t z 3> Tq. In the present paper we calculate the 
longitudinal interlayer magnetoresistance at Tq < At z < fiw c , 



which generalizes the result of Refs. [38|, l39|. This extends 
the applicability of the square-root dependence R zz oc \fB~ z 
to much wider region huj c > y/TQhu c > t z and gives linear 
Rzz (B z ) dependence in the interval ftiu c > At z 3> \ZTqHuj c . 

At 2t z > Tq one can use the 3D electron dispersion in Eq. 
([T|). Then the interlayer electron conductivity can be eval- 
uated at finite temperature using the 3D Kubo formula, 42 1 
which gives 

Ozz = J de [-n' F (s)} a zz (e), (3) 

where the derivative of the Fermi distribution function (43| 

n' F (e) = -l/{4Tcosh 2 [(s - n)/2T]}, (4) 
and the zero-temperature conductivity at energy e is 

= ^J2 y2 M [2lmG R (m,e)} 2 , (5) 



2 



where e is the electron charge, v z = de 3 o/dk z = 
2t z dsm(k z d) /H is the electron velocity, the sum over the 
electron quantum numbers m = {n, k y , k z } (excluding spin) 
is taken in the unit volume, and the retarded electron Green's 
function 



G R = [e-e 3D (™)-Z R (e,m)]~ 



(0) 



Both e3£)(m) and the electron self-energy is independent of 
k yi and the summation over k y in Eq. ([5]) gives the fac- 
tor equal to the LL degeneracy of one conducting layer per 
spin state Nll = l/2nl 2 { = e£? 2 /27rfi,c.[44| The interlayer 
conductivity is now given by a sum over LLs, 



o zz {e) 



&n (e) 



where the contribution to a zz from the n-th LL is 



e 2 hN LL 



2tt 



/ ^v 2 z (k z )[2lmG R (n,k z ,e)Y 



(7) 



(8) 



The electron self-energy E"(e, m) = E^(e) depends only 
on energy e and, possibly, on LL number n, being indepen- 
dent of k z (see Appendix). Substituting Eq. ^ and per- 
forming integration over k z , one obtains 



2nt 2 z |ImE«( £ ; 



-Re 



At 2 z - (Aef +iAe |lmE£(e)| 



4t*-(Ae-i\bn2X(e)\y 



> (9) 



where Ae = e — e 2 _o (n) — ReE^ (e), and cto denotes the in- 
terlayer conductivity without magnetic field: 



CO 



e 2 PF (v 2 z ) t = 2e 2 N LL t 2 z d/h 2 uj c r , 



(10) 



Pf = 2N[ J [ J /hijj c d is the 3D DoS at the Fermi level in the 
absence of magnetic field per two spin components, tq = 
h/2T and (v 2 ) = 2t 2 z d 2 /h 2 . 

To calculate the electron self-energy E^(e) we start from 
the standard model of 3D strongly anisotropic metals. The 
Hamiltonian consists of two terms: H = Hq + Hj. The 
first term H = Yl m ( TO ) c m c m describes 3D nonin- 
teracting electrons in magnetic field with anisotropic dis- 
persion given by Eqs. (HJ),©- The second term Hi = 
Si Vi ( r ) ^ + ( r ) ^ ( r ) describes the electron interaction with 
impurities. The impurities are taken to be point-like with 
the potential V$ (r) = US 3 (r — n) and randomly distributed 
with volume concentration n,. Without magnetic field the 
broadening Tq =ImE of electron levels due to the scat- 
tering by impurities in the Born approximation is Fo ~ 
irniU 2 pF, where pf is the density of electron states (DoS) 
at the Fermi level. In Q2D metals in strong magnetic field, 
fruj c > 4t z , r , the electron self-energy E^ depends on the en- 
ergy deviation Ae from the n-th LL. In the "non-crossing" 
approximation, 4|| 4(| schematically shown in Fig. [6j the 
electron Green's function, averaged over impurity configura- 
tions, has the form (see Appendix) 



G(ri,r 2 ,e)= ^ 



e - e 3D (n, k z ) - E n (e) ' 



(11) 



where k k ^ (r) are the electron wave functions in mag- 
netic field without impurities and E„ (e) is the electron self 
energy, averaged over impurity positions and given by the 
equation 



E n {e)=mU/[l~UG (e)] 



(12) 



Here G{e) = G(r,r,e) after the integration over k y and k z 
is given by 



N LL /d 



^ J(e-e 2D (n)-£„ (e)f~4t 2 z 



(13) 



The system of equations (|12|) and ([13]) allows to calculate 
the electron self-energy E n (e) numerically. 

For weak impurity potential UG (e) ~ UpF <C 1 the self- 
consistent Born approximation (SCBA) is applicable, and 
Eq. (jT2"j) reduces to 



E(e) w riiU + ni U 2 G (e) . 



(14) 



The SCBA is valid when the scattering potential of each im- 
purity is weak compared to Fermi energy, but the concentra- 
tion of impurities can be arbitrary, so that the broaden ing of 
electron levels Tq sw riiU 2 pf is also arbitrary. In Ref. [39( it 
was shown that one needs to apply at least SCBA to obtain 
qualitatively correct monotonic growth of interlayer longitu- 
dinal MR, which differs in SCBA and in the non-crossing 
approximations only by a factor close to unity. 

The similar system of SCBA equations for the electron 
Green's function in Q2D metals in magnetic field was written 
previously (see Eqs. (11) and (15)of Ref. [47j], Eqs. (11)- 
(13) in Ref. [H], Ec^ (29) in Ref. Eq. (17) in Ref. [3] 
or Eq. (6) in Ref. [50(). There the SCBA equations have not 
been solved, and the simple Born approximation was applied 



to calculate the interlayer conductivity. In Refs. [24j,|47J only 
the limit 2nt z S> huj c has been considered, while in Ref. [48| 
a large electron reservoir was introduced to damp the MOO 
and to make the simple Born approximation applicable. [5l| 
Consider the strong magnetic field limit, when hu> c > 
4t z ,To and the LLs do not overlap. Then one can consider 
only one LL in the solution of SCBA Eqs. flTJ) and ([13)1. 
which simplify to 



riiU 2 g LL /d 



(Ae - E,) 



4*5 



(15) 



where E* = E„(e) — riiU and Ae = e — 62d M — riiU = 
Ae+RcEf . At t z = we obtain the 2D result tea for the 
electron self-energy in SCBA. Note that the scattering en- 
ters Eq. (|15[) only in the combination riiU 2 gLL/ ! d = ToHuj c . 
Therefore, at t z 
in Eq. 



— > 0, T* = ^/TqTvjOc 
dT5J) , and naturally ImE» ~ ^/T^huic cx \[Bz accord- 
ing to Ref. [52j, which leads to the similar field dependence 
of background MR:^1 Hl| R zz w R zz0 |ImE (ji, B)\ /T oc 
\fB~ z . Eq. (fl"5"T) rewrites as an algebraic equation of the 
fourth power: 



is the only energy scale 

b: 



(Ae 



4tt 



(TofkjJcY 



(16) 
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FIG. 1: Imaginary part ImE of electron self energy as function 
of energy deviation Ae from the center of LL at t z /T* = (blue 
solid lme),t z /r» = 0.5 (green dashed line), and t z /T, = 1.0 (red 
dotted line). 
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FIG. 2: Real part of electron self energy ReE (Ae) at t z /T, 
0,0.5,1.0. 



Among four solutions of this equation, only one satisfies the 
physical requirement E* — > at Ae — > ±oo. Its solution 
is shown in Figs. 1 1 121 for three different values of 
i z /r* = 0,0.5, 1.0. One can see from Fig. 1 that ImE* 
in the finite energy interval of the width ~ 4r* at T* > 2t z 
and of the width ~ At z at T* < 2t z . ReE n has cusps at the 
boundaries of this energy interval. 

Now we substitute these solutions into Eq. ([9]) . The result 
for monotonic part of intcrlayer MR R zz = a zz is shown in 
Fig. [31 where a zz is the conductivity averaged over period 
fojj c of MQO. From this figure we see the crossover from 
linear to square-root field dependence of interlayer back- 
ground MR R zz (B z ). The interval of the linear MR is 
4t z < hjj c <C {4:t z ) 2 /Tq and increases with the increase of 
t z /To. The crossover from linear to square- root dependence 
of MR is a general feature of quasi-2D metals and already has 
been observed in a number of experiments (see, e.g. Refs. 
'MM 



The square-root dependence R zz oc B z is obtained 
at hui c > 2^/TqHujc > At z , which is much wider than the ap- 
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FIG. 3: Average interlayer magnetoresistance R zz (B z ) = l/a zz 
calculated from Eqs. and (|16l) as function of magnetic field 
at four different values of t z /To = 1 (solid green line), t z /To = 2 
(dashed blue line), t z /Fo = 3 (dotted red line) and tz/To = 4 
(dash-dotted black line). In the interval 4t z < Hui c <C (4t z ) 2 /To 
the interlayer MR has linear field dependence. 
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FIG. 4: Average interlayer magnetoresistance R zz (B z ) = l/a zz 
calculated from Eqs. © and (|16p as function of the square root 
of magnetic field at four different values of t z /To — 1 (solid green 
line), t z /To — 2 (dashed blue line), t z /To = 3 (dotted red line) 
and t z /To = 4 (dash-dotted black line). The linear dependence 
R zz (VB Z ) is obtained at huj c > 2\/hw c To > t z . Only in the 
interval 4t z < hw c (4f z ) 2 /To the interlayer MR has linear field 
dependence. 



plicability region hcu c 3> To ^> t z of the calculation in Refs. 

This result on the field dependence of interlayer MR can 
be analytically obtained as follows. Eq. (p~6|) simplifies to 
a quadratic equation at Ae = 0, i.e. at the center of LL, 

and has two solutions: E* = 2t 2 z ± ^ At\ + (T huj c ) 2 . The 
physical solution does not diverge at t z —> oo, has nonzero 



4 



imaginary part and, at t z = 0, agrees with the 2D limit 
described in Refs. 52j,|55[. All these criteria are satisfied for 
sign " — " , which at Ae = gives Re£* = and 



czz(e) 



llmE* 



4tj + (Tofuocf - 2tl 



(17) 



Substituting this to Eq. © one obtains the expression for 
the zero-temperature conductivity in the center of LL, i.e. 
at Ae = 0: 



(0) « 2a /ir. 



(18) 



In the limit Huj c > At z ^> \/Tohui c each LL gives an es- 
sential contribution to conductivity in the interval — 2t z < 
Ae < 2t z , as follows from Eq. (flB)) and can be seen from 
Fig. [5j where the conductivity as function of energy dis- 
tance from the nearest LL has a dome shape of the width 
~ 4t z . The conductivity, averaged over period fkv c of MQO, 
is then given by 



fhuj c /2 



-^o zz (Ae) « a zz (0) - — 



= a (2t„/hu c ) oc l/B 



7T U, 

4 huj c 



(19) 



This predicts a linear background magnetoresistance R zz = 
l/a zz oc B z in the interval At z < ftw c <C (4t z ) 2 /To of mag- 
netic field in quasi-2D strongly anisotropic compounds. The 
extra factor 7r/4 in Eq. (fT!?]) comes from the integration of a 
semicircle. Eq. (|T^1) also predicts stronger dependence of a zz 
on t z : a zz oc t\ unlike the usual dependence a zz ps cto oc t 2 ,. 

In the opposite limit fiu; c 3> Vmi^c ^ the "width 
of conducting band" from each LL is s» 4vmjHu^, and for 
background interlayer conductivity one obtains a zz /ao s» 
/BT, or 



2\/ro/Swc oc 1/ 



(0) « V , ^c/4r cx ^B z 



(20) 



in agreement with Refs. [38l - l40j and Fig. |4] 

The magnetic field B cr of the crossover from linear to 
square- root dependence of MR R zz (B z ) can be used to esti- 
mate the value of the interlayer transfer integral in the com- 
pound: t z ~ \/T\)hLj c , where uj c — eB cr /m*c corresponds to 
the crossover field. 

The obtained monotonic dependence of MR originates 
from MQO but survives at much higher temperature, when 
MQO are completely suppressed. According to Eqs. ©- 
((Sj, the temperature smearing of Fermi distribution func- 
tion does not influence the monotonic dependence of MR. 
Only at much higher temperature, when the electron scat- 
tering by phonons plays the major role in the relaxation 
of electron momentum, the above monotonic dependence of 
MR is weakened. Similarly, the long-range disorder, not in- 
cluded in the above calculation, produces the local variation 
of the Fermi energy along the sample on the scale greater 
than magnetic length, which damps the MQO but leaves the 
monotonic part of MR almost unchanged. Therefore, even if 
MQO are not seen in a compound because of strong disorder 
or high temperature, the background longitudinal interlayer 
MR can be observed and provides useful information about 
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FIG. 5: Interlayer conductivity <j zz (e) calculated from Eqs. ((9| 
and (|16p in the strong-filed limit hw c /To = 40 as function of 
energy counted from the nearest LL for four different values of 
tz/To = 0.1 (solid blue line), t z /To = 2 (dashed green line), 
tz/To = 5 (dotted red line), and t z /Fo = 10 (dash-dotted ma- 
genta line). In the limit 2t 2 z ^> TqFiw c the interlayer conductivity 
has approximately a semicircle shape of width 2t z . In the oppo- 
site limit Vofiujc 3> t 2 z the "conducting band" of each LL has the 
width ~ y/Yohw c . 



the electronic structure of this compound. This feature of 
background MR is similar to that of slow oscillations in Q2D 
compounds at t z > hujr.\22q 



To summarize, the calculation of longitudinal interlayer 
magnetoresistance in Q2D metals is performed in the strong- 
field limit Hu! c > At z . It predicts a linear background MR 
Rzz {B z ) — 1/ o zz oc B z in the interval 4i 2 < fku c <C 
(4t z ) 2 /Tq of magnetic field, where a zz has also unusual 
dependence on interlayer transfer integral: a zz oc t\. At 
stronger field or at smaller t z < TqHui c , the usual dependence 
a zz oc t 2 z is recovered, and the MR transforms to square-root 
dependence R zz (B z ) oc \fB~ z ^ a s was r ecently obtained in the 
limit t z < Iq < huj c in Refs. 



38 4 



The present calcula- 
tion generalizes these results to the region t z > To. The mag- 
netic field of the crossover from linear to square-root depen- 
dence of R zz (B z ) allows to estimate the interlayer transfer 
integral t z ~ yJTohuj c from experimental data. The obtained 
longitudinal MR R zz (B z ) explains numerous experiments on 
interlayer MR in strongly anisotropic quasi-2D compounds. 
The measurement of the monotonic part of longitudinal in- 
terlayer MR is much easier than the measurement of MQO 
or AMRO, because finite temperature and long-range crystal 
imperfections do not affect R zz (B z ) up to much higher tem- 
peratures or disorder. Therefore, the experimental study of 
longitudinal interlayer background magnetoresistance is pro- 
posed as a simple additional tool to investigate the electronic 
structure of strongly anisotropic quasi-two-dimensional com- 
pounds in a wide range of parameters. 

The work was supported by LEA ENS-Landau exchange 
program and by RFBR. 



APPENDIX 

Now we proof by the method of mathematical induc- 
tion, that in the "non-crossing" approximation the electron 
Green's function, averaged over impurity configurations, has 
the form of Eq. (fTTj) : 

r*ti \ ^n,k v ,k z { r l)^n,k v ,k z ( r2 ) . . 

G(ri,r 3 ,e)= > - -, . v -, , , 21 

nkk £ ~ e 3D{n,k z ) -T.nie) 

with E n (e) being independent of k v and k z . 

In the Born approximation with one impurity the electron 
self energy E^ (e) does not depend on k y and k z , but depends 
on the difference Ae„ = e — e 2 D (n) : 



a 



E r f ( £ )«C/ 2 G (e) = X; 



Nll/ d 



^{e- e 2D (n) - i 0) 2 - 4^ ' 



and its real part is an odd function of Ae„, ReE^ (Ae„) = 
— ReE^ (— Ae„), while its imaginary part is an even function 
of Ae„. Hence, in the Born approximation E„ (s, n,k yi k z ) — 
E„(Ae„) = E„(e). 

Assume that our condition is satisfied when each irredu- 
cable self-energy part contains no more than j impurities. 
Then the Green's function in coinciding points depends only 
on energy, Gj (e, r, r) = G (e): 



G 3 {e,r,r) = ]T 



E 

n,k z 



K* ky ,kS r )K.k y .kS r ) 

£ - £30 (n, k z ) - E n ,j (e) 

N LL /d 
£ - e 3 D (n, k z ) - T,„ tj (e) ' 



(22) 



When we include one more impurity to the irreducable elec- 
tron self-energy and average over its position, it remains in- 
dependent of k y , k z and depends only on Ae n : 

E„ J+1 = Urn/ [1 - UG 3 (Ae„)] . 

The Green's function Gj+i {e,r\,T2) also keeps the form of 
Eq. (fTTj) . To show this, let us find the function 

G i+ i (e,ri,r 2 ) = G (£,ri,r 2 ) 
^UG (ri,r) G (r,r 2 ) 



1 - C/Gj (e,r,r) 
Jd 3 rJ cPr'UGp (n,r) G (r, r>) G (r^) 
[1 - C/Gj (e, r, r)] [1 - UG 3 (e, r', r')] 

Go(£ ' ri ' r2)+ l-^G j(£ ) X 

K, ky ,kM) 



£ - (n,k z ) - E 0: „ (e) 
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FIG. 6: The set of diagrams for the irreducible self-energy, corre- 
sponding to the self-consistent one-site approximation. The dou- 
ble solid line symbolizes the exact electron Green's function. 



The integration over d 3 r in the last line gives 
S (n — n') S {k y — k' y ^ S [k z — k' z ), and the function 
Gj+i (e,ri,r 2 ) becomes 



G (e,ri,r 2 ) 

E 



Urii 



1 - UGj (e) 
K,k v ,kM)^° n * ky ^(r 2 ) 

k v ,k z t e ~ £ 3D (n, k z ) - E ,„ {e)\ 



+ .. 



i=0 



E 



• - e 3D (n, fc z ) - E ,„ (e) 

[1 - UGj (e)] 
e - e 3 D (n, fc z ) - E 0i „ (e) 

K,k y ,kM)^nX,kSr 2 ) 
■ - e 3D (n, k z ) - E i+ i jn (e) ' 



(23) 



where the new self-energy part 

E i+ i, B (e) = E ,„ (e) + Un t / [1 - UGj (e)] 

depends only on energy e and on LL number n. The Green's 
function in Eq. (1231 is again of the form of Eq. (ITT|) . 

The difference of impurity scattering to the same and to 
other LLs may give additional dependence of E„ (e) on LL 
number n, which requires further study. If E n) o (e) depends 
only on Ae„, then Gj (e) entering above formulas also de- 
pends only on Ae„, and E n j+1 (e) is also a function of Ae n 
only. If E nj o (e) depends only on e but not on n, and if the 
impurity scattering to the same and to other LLs has the 
same matrix element, then Gj (s) also depends only on e, 



and so does E 



n,j+l 



(e). 
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